Draft version February 2, 2008 
Preprint typeset using IiT^C style emulateapj 



EVIDENCE FOR A LINK BETWEEN FE Ka EMISSION LINE STRENGTH AND QPO PHASE IN A BLACK HOLE 

J. M. Miller 1,2 and J. Homan 3 

Subject headings: Black hole physics - relativity - stars: binaries - physical data and processes: accretion disks 

Draft version February 2, 2008 

ABSTRACT 

In X-ray binaries, the frequencies revealed in X-ray quasi-periodic oscillations (QPOs) are often interpreted 
as characteristic frequencies in the inner accretion disk, though the exact oscillation mechanism is unknown at 
present. Broadened Fe Ka lines are also excellent diagnostics of the inner accretion flow, if their broadening is 
indeed due to inner disk reflection. Herein, we present two cases where the flux and equivalent width of the Fe Ka 
emission lines in spectra of the Galactic black hole GRS 1915+105 vary with the phase of strong 1 Hz and 2 Hz 
QPOs in the X-ray flux. These results provide strong evidence that both QPOs and the Fe Ka lines originate in the 
inner disk. If the 1 Hz QPO is only a Keplerian orbital frequency, the QPO comes from a distance of 84 ± 26 Rsdm. 
from the black hole; the 2 Hz QPO implies a radius of 50 ± 15 Rschw At these radii, relativistic shaping of a disk 
line is inevitable. Moreover, the link holds in radio-bright and radio-faint phases, signaling that in systems like 
GRS 1915+105, the Fe Ka line is a disk line and not a jet line as per SS 433. A particularly interesting possibility 
is that a stable warp in the inner disk, e.g. due to Lense-Thirring precession, may produce the observed QPOs and 
line modulations. More broadly, the Fe K-QPO link provides an unprecedented mechanism for revealing the inner 
accretion flow and relativistic regime in accreting systems, in that it gives two measures of radius: for a given disk 
QPO model, the frequency translates into a specific radius, and relativistic line models yield radii directly. 



1. INTRODUCTION 

In accretion-powered sources, a great deal of theoretical and 
observational effort is devoted to the study of the inner accre- 
tion flow, and to using the accretion flow itself as a tool to study 
the central compact object. High frequency QPOs in Galactic 
black hole systems may be related to the Keplerian frequency 
at the innermost stable circular orbit (for zero-spin black holes, 
v = 220 Hz 10 Mq/M bh ; e.g., Strohmayer 2001, Miller et al. 
2001, Remillard et al. 2002, Homan et al. 2004), but they are 
only observed in a handfull of sources. QPOs at lower frequen- 
cies, typically 1-10 Hz are much more common. Although 
their origin is unclear, interpreting them as a Keplerian fre- 
quency puts them at <;100 Rschw.- If these QPOs are a preces- 
sion frequency at the inner disk (e.g., due to frame dragging, 
see Markovic & Lamb 1998), they probe regions much closer 
to the black hole. Broad, asymmetric Fe Ka emission lines 
are found in the spectra of both supermassive black holes and 
stellar-mass Galactic black holes (Reynolds & Nowak 2003). 
These lines appear to be imprinted with the extreme Doppler 
shifts and gravitational red-shifts expected near to a black hole, 
but the case is not entirely shut. 

A natural evolution of the means by which inner accretion 
flows are studied is to bring spectral and timing constraints to 
bear jointly. It has been shown that the disk reflection spectrum 
(the Fe Ka line is the most prominent part of the reflected spec- 
trum; see George & Fabian 1991) in Galactic black holes varies 
with broad-band frequency (Gilfanov, Churazov, & Revnivtsev 
2000). This is strong evidence that timing and spectral proper- 
ties are intimately related, and supports models wherein Fe Ka 
emission lines arise in the inner disk. 

In this work, we take spectral and timing connections a step 



further. We report the discovery of a link between discrete 
timing features (QPOs) and spectral features (Fe Ka emission 
lines): Fe Ka emission line strength varies with QPO phase 
in radio-bright and radio-faint states of the "microquasar" 
GRS 1915+105. Our results provide the first indication that 
these diagnostics of the inner accretion flow are related. 

2. DATA SELECTION AND REDUCTION 

We chose to investigate GRS 1915+105 as strong QPOs (see, 
e.g., Morgan, Remillard, & Greiner 1997) and Fe Ka emis- 
sion lines (Muno et al. 2001, Martocchia et al. 2002, Miller et 
al. 2004a) have been reported in this source. The high mass 
of GRS 1915+105 (M BH = 14.4±4.4M Q , Harlaftis & Greiner 
2004) means that for a given QPO frequency, the QPOs may 
trace regions closer to the black hole in this source than in 
lower-mass black holes. GRS 1915+105 displays a wide va- 
riety of radio states. Steady radio emission is associated with 
compact jets (Fender 2004), and this fact can be exploited to ex- 
amine whether any relation between Fe Ka line flux and QPOs 
depends on the presence of a jet. 

To extract QPO phase-selected spectra, QPOs should be 
strong enough to dominate the X-ray variability and be di- 
rectly visible in the X-ray lightcurves. Morgan, Remillard, 
& Greiner (1997) found large amplitude, sinusoidal QPOs in 
GRS 1915+105 that could be seen in the naked lightcurve. 
The QPOs were found to be quasi-periodic largely due to 
random jumps in phase. We selected two observations of 
GRS 1915+105 available in the RXTE public archive: 20402- 
01-15-00 (obs. date: 9 Feb. 1997) and 20402-01-50-01 (obs. 
date: 16 Oct. 1997). The spectral and timing properties of these 



observations (including an Fe Ka emission line) 
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FlG. 1 . — The power density spectra of the radio-bright observation (black) 
and radio-faint observation (red) are shown above. The strongest QPOs have 
frequencies of 1.052(3) Hz and 2.260(2) Hz, and amplitudes of 14.7(3)% rms 
and 13.6(1)% rms, respectively. 

have previously been reported on by Muno et al. (2001). Im- 
portantly, Muno et al. (2001) also studied the radio properties 
of GRS 1915+105 during these observations; the first occurred 
in a "radio-faint" state (hereafter, "RF"; the source was not de- 
tected at 8.3 GHz), and the second in a radio-bright or radio- 
plateau state (74 ± 6 mJy at 8.3 GHz; hereafter "RB"). 

Tools within LHEASOFT version 5.3 and CIAO version 3.1 
were used to reduce and analyze these observations. Standard 
RXTE data screening and good time selections were applied 
(e.g., for SAA avoidance). We further screened the data to only 
include times during which all five PCUs were operated. The 
net total good time for the RF obs. (15-00) was 10.1 ksec, and 
the net total good time for RB obs. (50-01) was 4.5 ksec. All 
spectra and lightcurves were generated using the tool "saex- 
trct", and all response matrices were generated with the tool 
"pcarsp". As GRS 1915+105 is very bright and we are inter- 
ested in relative variations, background spectra were not sub- 
tracted. We made time-averaged broad-band PCA spectra by 
combining data from all layers of all PCUs taken in "Standard- 
2" mode (129 energy channels between 2 and 60 keV, taken 
every 16s). Spectra and lightcurves were also made from data 
taken in "B_8ms_16A_0-35_H" mode (16 energy channels up 
to channel 35 - roughly 13 keV - taken every 8ms). 

To accurately estimate relative flux differences in data ac- 
quired with the same detector, it is only important to under- 
stand the degree to which systematic errors may change over 
time. For this purpose, we reduced and analyzed Crab spectra 
obtained before and after each observation of GRS 1915+105 
(spectra obtained on 31 Jan. 1997 and 16 Feb. 1997, and 12 
Oct. 1997 and 27 Oct. 1997). Absorbed power-law fits to these 
spectra reveal that most channels below 10 keV (those most 
important for Fe Ka line studies) differ by 0.2% or less. To be 
conservative, we added 0.2% systematic errors to the spectra. 
This is less than the 0.5-1.0% systematic errors often added to 
PCA spectra when absolute fluxes are of interest. 

3. ANALYSIS AND RESULTS 

We made power-density spectra (PDS) of the X-ray flux from 
each observation using all available data, and fit the 
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FlG. 2. — Approximately 20 cycles of the QPO wave from randomly 
selected intervals in each observation are shown above. Clearly the QPOs 
strongly dominate the lightcurves. Solid lines mark the mean count rate, and 
dashed lines above and below mark the count rate levels used to select the high 
and low phases of the QPO waves. The radio-bright data are binned to 0.128 s, 
and the radio-faint data are binned to 0.064 s. 

0.01-100 Hz range with Lorentzian functions in the j/-max rep- 
resentation (see, e.g., Homan et al. 2004). Each PDS is strongly 
dominated by a single QPO (see Figure 1). The strongest QPO 
in the RF obs. has a frequency of 2.273(4) Hz, a (2-value of 
6.0(2), and an amplitude of 11.8(2)% rms (errors are lcr). The 
strongest QPO in the RB obs. has a frequency of 1.050(3) Hz, a 
Q-value of 7.9(4), and an amplitude of 1 1.1(2)% rms. Assum- 
ing ablack hole mass of Mbh= 14.4±4.4M Q , if these frequen- 
cies are merely Keplerian orbital frequencies, they correspond 
to radii of 84 ± 26 Rsdm. ar, d 50 ± 15 Rsdiw. respectively. 

Lightcurves from the "B_8ms_16A_0-35_H" mode data 
from each observation were analyzed, and the total mean count 
rate and mean count rate in 100s intervals were calculated. The 
mean rate in the RF obs. is 3532 c/s, with a range of 3298-3727 
c/s when measured in 100s intervals. The mean rate in the RB 
obs. is 6602 c/s, with a range of 6496-6793 c/s. 

As the QPOs in each observation dominate the X-ray vari- 
ability (see Figure 2), and because the mean count rate in each 
observation is very steady, we applied simple count-rate selec- 
tions to isolate the maxima and minima of the QPO waves. For 
the RF obs. selecting periods when the flux was 15% above 
or below the mean effectively isolated the maxima and minima. 
For the RB obs. there is less noise apart from the QPO, and 
periods when the flux was 10% above or below the mean were 
found to effectively isolate maxima and minima. The CIAO 
tool "dmgti" was used to generate additional good time files of 
the maxima and minima. These files were than applied within 
"saextrct" to produce spectra of the maxima and minima. 

All spectral analysis was done using XSPEC version 11.3 
(Arnaud 1996) and IDL version 5.4. Analysis of the time- 
averaged standard-2 spectra was performed in the 3.0-20.0 keV 
band. The PCA is calibrated well in this energy range, and this 
range is not much greater than the range of the "B_8ms_16A_0- 
35_H" mode data. Fits to the spectra with a number of common 
models revealed the Fe Ka emission line previously noted by 
Muno et al. (2001). In each case, the addition of a simple 
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FlG. 3. — Fits to the high and low phase spectra from the radio-faint ob- 
servation of GRS 1915+105 are shown above. The spectra were fit with the 
model described in Table 1, using the fixed line centroid and width (which 
gives a smaller difference in the line fluxes than if these parameters float). The 
line flux was set to zero on the ratio plot to show the difference in Fe Ka line 
intensity between high and low QPO phases. 

Gaussian component to model the line was significant at more 
than the 8er level of confidence. We find that disk reflection 
models provide the best overall fit to the broad-band spectrum 
and do not require a disk that is prominent in the RXTE band- 
pass. Our model does not include a disk, but this only reflects 
the fact that a cool disk with moderate flux is not easily seen in 
the RXTE bandpass, especially when Nh is high. 

We fit the time-averaged, high phase, and low phase spectra 
with absorbed cut-off power-law (for the radio-bright phase) 
or broken power-law (for the radio-faint phase) models, mod- 
ified by the addition of a Gaussian component (to model the 
Fe Ka emission line) and smeared edge component ("smedge" 
in XSPEC) to mimic a reflected continuum. Trudolyubov 
(2001) also found that these differing spectral models were re- 
quired in RF and RB states. The "phabs" model was used to 
account for absorption in the neutral ISM. The "smedge" was 
not statistically required in all fits, and its parameters could not 
be well constrained. However, the smedge effectively mimics 
a disk reflection spectrum at energies above the Fe Ka emis- 
sion line, so we fixed the model parameters to modest values 
consistent with the data (E = 8.5 keV, r = 0.3, width= 10 keV). 

The parameters obtained from the spectral fits are shown in 
Table 1 . In all cases, the line fluxes and equivalent widths are 
significantly different in the high and low phase spectra (see 
Fig. 3). The confidence intervals do not overlap even when the 
ltr error bars are multiplied by 5. There is evidence that the 
Fe Ka emission line flux is broader and slightly shifted to the 
red in the high phase relative to the low phase. 

To demonstrate that the line flux varies in a model- 
independent manner, we subtracted the low count rate spectra 
from the high count rate spectra, and divided the result by the 
mean count rate for the entirety of each observation. Using 
the energy to channel bounds found in the response matrix for 
each observation, the bins in the difference spectra were con- 
verted to energy bounds; however, the detector response was 
not removed and no spectral fitting was performed. In both 
the radio-bright and radio-faint observations, the most promi- 
nent difference between the high and low phases is the Fe Ka 
line (see Fig. 4). For a "control" observation, we extracted 
high and low phase spectra (10% above/below the mean) from 
a 1.2 ksec slice of an observation of the Galactic black hole 
H 1743-322 obtained on 28 May 2003. Like the observations 
of GRS 1915+105, H 1743-322 was highly variable on 
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FlG. 4. — The plot above shows the relative intensity of the radio-bright 
and radio-faint spectra of GRS 1915+105 (and a "control" observation of 
H 1743-322) in high and low count rate phases. The low count rate spectra 
were subtracted from the high count rate spectra, and the difference was di- 
vided by the mean count rate for the entire observation. Error bars are plotted 
but are very small. Channel bins were converted to energy space but the detec- 
tor response has not been removed (so the Xe L edge complex between 4-6 keV 
is seen). Unlike the H 1743-322 difference spectrum, the Fe Ka line bins are 
dominant in the GRS 1915+105 difference spectra. This model-independent 
analysis confirms that the Fe Kct line intensity varies with QPO phase. 

timescales of seconds, but had a steady mean count rate on 
100 s scales in this window. Spectroscopy of H 1743-322 in 
outburst did not reveal an Fe Ka emission line (Miller et al. 
2004b), and the Fe Ka line bins are not globally dominant in 
the difference spectrum shown in Fig. 4. 

4. DISCUSSION AND CONCLUSIONS 

The main results of this work may be summarized as follows: 

• Fe Ka line flux (and possibly FWHM and centroid energy) 
depends on QPO phase in 1-2 Hz QPOs in GRS 1915+105, 
linking two independent diagnostics of the inner accretion flow. 

• The QPO frequencies at which the Fe Ka line flux is modu- 
lated correspond to radii of 84 ± 26 Rs c hw. an d 50 ± 15 Rsdiw., if 
they are merely Keplerian orbital frequencies. This link likely 
ties Fe Ka lines in Galactic black holes to the inner disk and 
further supports evidence that the observed line broadening is 
due to dynamics in the disk (Doppler shifts, gravitational red- 
shifts) close to the black hole. 

• The presence of Fe Ka lines, QPOs, and the link between 
them does not depend on jet activity (taking radio flux as a jet 
indicator). Although the base of a failed jet may be a source of 
hard X-rays which irradiate the disk, Fe Ka lines and QPOs do 
not arise in an extended jet. 

Blobs in the inner accretion disk might create QPOs and 
could give a quasi-periodically changing reflector area. Sim- 
ilarly, if the height of hard X-ray emitters (flares, or the base 
of a jet) changes quasi-periodically, QPOs and Fe Ka line 
variations might be expected. However, while both explana- 
tions are possible, such explanations seem inconsistent with 
the strength and constancy of the QPOs found in these obser- 
vations of GRS 1915+105, and with the fact that the Fe K- 
QPO connection holds over a factor of two in flux and in both 
radio-bright and radio-faint phases. Moreover, associating low- 
frequency QPOs with Keplerian frequencies is problematic for 
two reasons: (1) X-ray production in accretion disks should be 
centrally concentrated, with little emission from r > 100 Rs c hw., 
and (2) QPOs are stronger at higher energies, suggesting they 
are also centrally produced (see, e.g., Homan et al. 2001). A 
stable warp at the inner disk, e.g. due to Lense-Thirring preces- 



4 The Iron Line - QPO Connection 

Table 1 

Fits to QPO Phase Spectra 
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NOTE. — "RF" refers to the radio-faint observation (fit with a broken power-law), and "RB" to the radio-bright observation (fit with a cut-off power-law). A 
"smedge" with E=8.5 keV, r = 0.3, and width=10 keV was included in all fits. Fits to "avg" spectra are fits to standard-2 data in the 3-20 keV band. All other fits are 
to lower resolution "B_8ms_16A_0-35_H" data in the 3-13 keV band, t Denotes fits in which the Fe Ka line centroid and FWHM were fixed to the values measured 
in the average spectrum. Where errors are not given, parameters from the standard-2 fits were fixed. Systematic errors of 0.2% were added to all spectra prior to fitting 
to estimate the drift in detector response. Fits to the phase-selected spectra with line centroid energy and FWHM fixed to the average value, and with parameters free, 
are both reported for comparison. All errors are 1<t errors to allow confidences to be inferred directly. Formally acceptable fits are obtained with systematic errors set 
to 0.6%, which is typical for spectroscopy wherein absolute flux measurements are important. 



sion, could cause QPOs and would also serve to create a quasi- 
periodically changing reflector area, modulating the strength of 
Fe Ka line emission. It has been suggested that Lense-Thirring 
precession may be the cause of low-frequency QPOs in accret- 
ing systems (e.g., Markovic & Lamb 1998), and the Fe K-QPO 
connection may be evidence for this General Relativistic effect. 

Whether or not Lense-Thirring precession can account for 
the Fe K-QPO connection, the Fe K-QPO connection has the 
potential to reveal the innermost regime in accreting systems. 
For a given disk-driven QPO model, the Fe K-QPO connection 
gives two measures of radius: frequencies imply specific radii, 
and relativistic disk line models constrain radii directly (Fabian 
et al. 1989, Laor 1991). The connection over-constrains the 
system, if QPOs are disk frequencies. 

Work of this kind in Galactic black holes may have impor- 
tant implications for Fe Ka line reverberation mapping in AGN 
with planned missions such as Constellation-X and XEUS (e.g., 



Young & Reynolds 2000), and observations of disk hot spots 
orbiting close to the black holes in AGN (e.g., Iwasawa, Mini- 
utti, & Fabian 2004). Long RXTE exposures with carefully 
chosen PCA modes are needed to obtain additional data for 
new studies. The large effective areas of XMM-Newton and 
Astro-E2 may support such studies as well, and may even allow 
studies of the line shape with QPO phase. A new X-ray mission 
with an area of 2- 10m 2 and CCD spectral resolution (or better) 
may fully exploit the Fe K-QPO connection. 
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